Anisotropic super-spin at the end of a carbon nanotube 
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Interaction-induced magnetism at the ends of carbon nanotubes is studied theoretically, with a 
special focus on magnetic anisotropies. Spin-orbit coupling, generally weak in ordinary graphene, 
is strongly enhanced in nanotubes. In combination with Coulomb interactions, this enhanced spin- 
orbit coupling gives rise to a super-spin at the ends of carbon nanotubes with an XY anisotropy on 
the order of 10 mK. Furthermore, it is shown that this anisotropy can be enhanced by more than 
one order of magnitude via a partial suppression of the super-spin. 

PACS numbers: 73.63.Fg, 75.70.Tj, 75.30.Gw, 75.50.Xx 



Introduction. Carbon-based material systems are 
promising candidates for applications in future (quan- 
tum) information technologies and attract much atten- 
tion. One reason for their great potential is that elec- 
trons on carbon honeycomb lattices behave unusually in 
many respects. For instance, instead of being quasiparti- 
cles resembling free electrons, the low-energy excitations 
in graphene and in some carbon nanotubes (CNTs) are 
described by the Dirac theory [l(. An important con- 
sequence of this Dirac nature is the efficient suppres- 
sion of Coulomb interaction effects, even though, due to 
the extreme electronic confinement to only one atomic 
layer, the energy scale corresponding to this interaction is 
rather large on two-dimensional carbon honeycomb lat- 
tices Q. As a consequence, non- interacting theories of 
graphene or CNTs are very successful as long as the hon- 
eycomb lattice is intact; in order to drive a phase transi- 
tion in bulk graphene, even larger interaction strengths 
are required |3(. 

Where the honeycomb lattice is disturbed, however, 
the Coulomb interaction has a chance to manifest itself. 
Examples of this mechanism are the vacancy-induced 
magnetism in graphite [H or edge magnetism, which ap- 
pears at zigzag edges of graphene flakes In these 
examples, disturbances of the honeycomb lattice (vacan- 
cies and edges, respectively) allow interactions to drive a 
magnetic transition. In the present work, a similar mag- 
netic transition at the ends of CNTs @ is studied: the 
Coulomb interaction aligns the spins of all electrons in 
certain states, localized at the CNT end, and thus gives 
rise to a super-spin, composed of many individual elec- 
tron spins. This CNT end magnetism is the pendant of 
graphene's edge magnetism 

Usually, the spin-orbit interaction (SOI) is assumed 
to be negligible for edge magnetism. This is because 
the SOI energy scale in graphene (28 /ieV Q) is much 
smaller than the electron-electron interaction (~ eV Q). 
Moreover, the rigidity of the super-spin suppresses SOI 
effects even further so that their effective energy scale 
is in the jjK regime, and this renders the SOI experi- 
mentally irrelevant. However, it is known that the SOI 
can be enhanced: at graphene/graphane interfaces, for 



instance, spin-orbit effects have been shown to be am- 
plified by two orders of magnitude ||. Also the surface 
curvature of CNTs gives rise to an enhanced SOI (loj . 

In this Letter, the spin-orbit anisotropy of the super- 
spin at CNT ends is calculated on the basis of a micro- 
scopic model. It is shown that the curvature-induced 
enhancement of the SOI lifts the anisotropy up to an 
experimentally accessible regime. Moreover, by a com- 
bination with the mechanism for tuning the strength of 
edge magnetism [ll[ , which translates to a tuning of the 
size of the super-spin at CNT ends, the anisotropy can 
be further increased. 




FIG. 1. (Color online) Left: the end of a 15-CNT, i.e., a CNT 
with N c = 15 unit cells in circumferential direction. Right: 
the local coordinate system {f,t,z} (z points into the paper 
sheet) and the orientation of the p orbitals p r and pt at the 
surface of the CNT (dashed line). 

Model. This study starts from an interacting model 
for the 7T and a orbitals of the CNT's carbon atoms. The 
Hamiltonian consists of three parts 



H = H hop + H 



so 



H 



int 



(i) 



The hopping between different orbitals of nearest- 
neighbor atoms is described by H hop = $ ,c\ Cjn'e + 



where the fermionic operator annihilates 



an electron with spin a in orbital fj, at site i of a hon- 
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eycomb lattice rolled up to a CNT. Here and hence- 
forth, repeated indices are to be summed over. t l £ , is 
the amplitude for an electron hopping from orbital p! 
at site j to orbital p at site i. All second shell or- 
bitals p = s,p r ,pt,p z of the carbon atoms are taken 
into account. The p orbitals are defined in a local co- 
ordinate system {r,t,z} (see Fig. [I]). e s = —8.9 eV 
is the energy of the s orbital. The CNT curvature 
in a highly nontrivial manner, de 
la and ~ 



is encoded in t*?„. 
scribed in detail in Refs 
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The SOI reads 
with Ago = 6 meV 
is the Levi-Civita tensor and s M are the spin 
Pauli matrices (with eigenvalues ±1), defined in the lo- 
cal coordinate system {r, t,z}. The contributions from 
the d orbitals to the SOI are negligible because of the 
curvature- induced tto~ hybridization [7| . The Coulomb in- 
teraction is modeled by a Hubbard Hamiltonian for the 
7r orbitals only, i.e., H lnt = V ' c\ Pi ^Ci Pr ^c\ pri c iPr ^, with 
U ~ 9.3 eV 0. The Coulomb interaction for the a elec- 
trons is assumed to be included on an effective level in 
the hopping parameters. 

Edge state model. The CNT considered in this work 
is half-infinite, i.e., it has only one end at which the lat- 
tice is zigzag-terminated (see Fig. [1]). Such a termina- 
tion supports edge states [l^. For describing the mag- 
netic state at this CNT end, it is sufficient to work in 
a restricted Fock space, spanned by these edge states 
only [ll|. Therefore, fermionic edge state operators 
efco- = n)]*c-in<! are introduced, where ip k (i,p) is 

the edge state wave function with momentum k in cir- 
cumferential direction (along t). ip k (i,n) is determined 
by the numerical diagonalization of the hopping Hamil- 
tonian H hop . It is derived mainly from the tt orbital 
p r . However, due to curvature- induced 7rer hybridiza- 
tion, it has also contributions from the a bands; these 
are most important for the SOI. Simple approximations 
for jjl) takin g in to account only the tt band can be 
found in Refs. la, LLjJ, and UM These approximations are 



sufficient for calculating the projection of the Hubbard 
Hamiltonian and of the kinetic energy of the edge states 
(see below), but they are insufficient for the SOI. 

Essentially, there are three terms describing the 
physics of the edge states. The kinetic energy of the 
edge states [HI is described by 



H. 



kin 



k e \a e kcr-> 



(2) 



ka 



with Afk = y2 cos(fc — 7r) — 1. The k sum is restricted 
to ^ < k < ^- since the edge states exist only for this 
range of momenta in circumferential direction |16j . A is 
a phenomenological parameter which accounts for a va- 
riety of electron-hole symmetry (EHS) breaking terms in 
the Hamiltonian, such as next-nearest neig hbor hopping 
or local perturbations at the edge [ll], Elf- With EHS, 
the edge states have exactly zero energy. Breaking EHS 



gives rise to a non-zero bandwidth for the edge states; in- 
terestingly, all these EHS breakings can be subsumed in 
one parameter A. Note that A is experimentally tunable 
in special geometries [ll| . 

The projection of the honeycomb lattice Hubbard in- 
teraction H int onto the edge state subspace reads [ll 



H" 



U 



— 2^ T(fc, k', q) : el +q ^el^ qi e kli :, (3) 



"k,k' ,q 



where U ~ 9.3 eV [8(] is the Hubbard coupling con- 
stant on the 2D honeycomb lattice, N c is the number of 
unit cells in circumferential direction, and T(k,k',q) — 
A/fe+gAfeA/fe'-gA/fc'/C 1 _ u* k+q u k u* k ,_ q u kl ), with u k = 1 + 
e tk . Again, the k sum is restricted such that the momenta 
of all edge state operators are in the interval ] ^ , [. 

The interplay of H^ ln and H™ 1 has been extensively 
discussed in view of the tunable ferromagnetism at long 
graphene zigzag edges 11, The present work rather 
focuses on finite zigzag edges with periodic boundary con- 
ditions (the graphene edge is rolled up to a CNT end). 
Thus, k is discrete with a spacing Ak = 2tt/N c . For 
A <C U, a super-spin with a certain size S is found at 
the CNT end and, because Hf n and arc SU(2) in- 
variant, this super-spin obeys the full rotation symmetry, 
that is, the ground state is (2S + l)-fold degenerate. 

The effective SOI is derived by projecting H so onto the 
Hilbert space spanned by the edge states. Since is a 
single particle operator this projection is easily performed 



(4) 



k,k'a,a' 



where \tp k , o~) is the wave function of an edge state with 
momentum k and spin a. H so can be separated into a 
part H SO ' z proportional to s 2 , and a part H SO ' xy propor- 
tional to s x or s y . Because H so depends on the locally 
defined s r ' f ' z , the terms H so,xy have a trigonometric de- 
pendence on the azimuthal angle <p, while H SO ' z is inde- 
pendent of <p 17j. Thus, the projected H^°' z is diagonal 
in k, while iff°' xy is not. 

The most important part of the effective SOI for the 
super-spin anisotropy is the /c-diagonal i/f° :Z . Due to 
time-reversal invariance, ae z (k) = (ipf., a\ H so \tp k , a) 
must be odd around k = and k = tt, so that e z (k) can 
be expanded in sin(m/c), with integer m. It turns out to 
be sufficient to keep only the leading order, i.e., e z (k) = 
5 Z sin(fc). The parameter S z is determined by fitting e z (k) 
to the numerically evaluated (V>|)t| H IV'feit)- As can 
be seen from Fig. ® the sme form for e z (k) is an excellent 
approximation so that one may write 



H S e ° z = S Z J2 sin(k)ael a e ka . 



k.a 



(5) 



The effective spin-orbit coupling constant S z is fitted to 
the numerical results for N c — 5, 60 in order to extract 
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the N c dependence of S z 



1.6 1.1 

Sz /A so c — + w 



(6) 



The relative deviation of Eq. © from the numerical data 
is less than 10~ 4 . 
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FIG. 2. (Color online) Numerical results (black dots) and fit 
S z sin(fc) (red line) for the effective z SOI of edge states. In 
this figure, N c = 10. The deviation of the numerical e z from 
the sine form at k = ^S- is due to finite size effects in the 



numerical calculation, 
for N c = 10. 



The inset shows the function e xy (k) 



The functional form of the Hamiltonian H^°' xy is less 
obvious. First of all, due to the trigonometric depen- 
dencies on the azimuthal angle, H^° xy is not diagonal 
in k but couples momenta differing by Ak = 2tt/N c . In 
general, one may write 



^jSO,xy 



k 



J fc+Afc/2 



'(*)«£ 

e*v(k) = Re[(^ +Afe/2 ,t l# SO |^._ 



-is v )-e k _ Ak/2 + R.c. 

(7) 

Afe/2!^)] is a g am de- 



termined numerically and e k = (e^, e]L ). The result is 
shown in the inset of Fig. [5J Note that e xy (k) is signifi- 
cantly smaller than e z (k). Furthermore, since the super- 
spin which is generated by the Coulomb interaction has 
zero momentum in circumferential direction, the k non- 
diagonal iff °' xy is only of minor importance and will be 
neglected in the remainder of this work. 

Anisotropy of the super-spin. Usually, the kinetic en- 
ergy of the edge states, quantified by A, is small com- 
pared to the interaction energy scale U, so that the spins 
of all N ~ N c /3 electrons occupying the edge states are 
aligned (half-filling is assumed) and form a super-spin 
of length 5 = N/2 In the absence of SOI, this 

super-spin is rotationally invariant, i.e., the ground state 
is (25 + l)-fold degenerate. For the CNTs considered 
here, this degenerate subspace is protected by excitation 
energies of a few hundred meV, so that the magnetic 
properties of a CNT end are expected to be well described 
by an effective super-spin model living in this (25 + 1)- 
dimensional subspace. All operators in this super-spin 



model are given in terms of the (25 + l)-dimensional 
representation of the angular momentum operators S v , 
with v = x,y,z, for which £„(5^) 2 = 5(5+1). 
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FIG. 3. (Color online) Spin response of the super-spin state 
at the end of a N c — 17 CNT. The dots are the results of the 
exact diagonalization of the interacting edge state Hamilto- 
nian. The red line is the spin response calculated from the 
super-spin model Eq. (SJ. 

The effective super-spin Hamiltonian H ss is determined 
by requiring equal spin responses of the super-spin model 
H ss and the edge state model H e = Hf n + H*° '< z + H™ 1 to 
a Zeeman field H z = i.e., (5% = (S v ) e . 

Here, (-) ss denotes the ground state expectation value 
with respect to H ss +H z and (-) e the ground state expec- 
tation value with respect to H e + H z . In the edge state 
model, the spin operators entering the Zeeman Hamilto- 

,t 



nian are given by S v — | ^ e k 
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TABLE I. Maximal super-spin S, bare anisotropy constant 
D, enhanced anisotropy constant D* , and radius R of various 
zigzag CNTs with different jV c . 



Following Ref. [15|, the edge state model is diago- 
nalized exactly in order to calculate the spin response 
of the ground state. In the absence of SOI and Zee- 
man fields, the ground state is (25 + l)-fold degener- 
ate (S z = —5,..., 5). For Ago ^ 0, this degeneracy is 
slightly lifted. As the SOI proportional to S x ' v is ne- 
glected, S z is a a good quantum number also in the pres- 
ence of SOI, and the state(s) with the smallest \S Z \ forms 
the new ground state. For integer (half-integer) 5, the 
ground state consists of S z = (5 Z = ±-|) and is non- 
degenerate (two- fold degenerate). Thus, the SOI- induced 
anisotropy is of XY type, i.e., 



D(S z f 



(8) 
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with positive D, determined numerically. This is one of 
the main results of this work. 

For a 17-CNT one finds D = 0.11 fieV by fitting D to 
the position of the first step of the z-spin response. Fig. 
[3] shows that the effective super-spin model describes the 
exact spin response extremely well. For instance, the 
relative differences of the positions of the steps in the S z 
response is typically 10~ 6 . In the limit Ago <C U, the 
dependence of D on the SOI and Coulomb interaction 
strengths is D ~ Ag Q /[/. Table fl] presents an overview 
of anisotropy constants for different CNTs. 
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FIG. 4. (Color online) Increased anisotropy for the partially 
suppressed super-spin. For small |A| the size of the super-spin 
is S — 3 and for larger |A|, S = 1. The dashed vertical lines 
mark the transition between these regimes. The insets show 
schematically the configurations of the edge state electrons: 
arrows correspond to occupied electrons with a certain spin 
- light arrows are paired spins and bold arrows are unpaired 
spins; gray circles are empty states. The dashed parabolas 
symbolize the mean-field energies of the edge states [rD |. 



Enhanced anisotropy. The configuration in which all 
spins are aligned (for |A| <C U) corresponds to the satu- 
rated edge magnetism (SEM) regime [Iff - In this phase, 
the system is very rigid, i.e., the energy gap to the first 
excitation is large. This rigidity suppresses SOI effects. 
The weak edge magnetism (WEM) phase (for |A| ~ U) 
is less rigid and has smaller excitation gaps. It is thus 
expected that D is enhanced in the WEM phase. In the 
limit of infinitely long graphene edges, as it has been dis- 



cussed in Rcfs. |11| and |15|, it is possible to reduce the 



edge magnetization continuously. A CNT end, however, 
corresponds to a short edge and the size of the super- 
spin can be reduced only in steps of AS" = 2. Thus, the 
smallest CNT for which the super-spin can be reduced to 
S > 1 is the 17-CNT. The anisotropy enhancement will 
be discussed on the basis of this CNT, but Tab. U shows 
the enhanced anisotropics also for other CNTs with even 
electron numbers. The case of odd electron numbers is 
more complicated and beyond the scope of this work. 
The dependence of D on A is shown in Fig. [JJ For 



small |A|, D does not depend on A. Near the critical 
kinetic energy strength |A| ~ 1.3 eV, at which the re- 
duction of the length of the super spin from S = 3 to 
5 = 1 takes place, D increases by more than an order of 
magnitude. The maximal D is D* ~ 3.1/LteV. Compared 
to the anisotropy of the super-spin with the maximal size 
(SEM regime), this is an enhancement by a factor of 27. 

Conclusion. The SOI-induccd anisotropy of the super- 
spin at CNT ends has been calculated and it was shown 
that this anisotropy is of XY-type. Two mechanisms 
have been discussed which enhance this anisotropy. One 
is based on the well-known curvature-induced SOI en- 
hancement [l(| and the other is based on the tunability of 
edge magnetism 



11 



The first effect lifts the anisotropy 
of small CNTs to experimentally accessible regimes. The 
second effect can be used to increase the anisotropy also 
in larger tubes. These effects are independent. Thus, it 
is expected that the latter effect is also relevant in other 
SOI-cnhanced systems supporting magnetic edge states, 
in particular at graphene/graphane interfaces 
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